Optical analogue of population trapping in the continuum: classical and quantum 

interference effects 
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A quantum theory of light propagation in two optical channel waveguides tunnelling-coupled to a 
common continuum of modes (such as those of a slab waveguide) is presented, and classical and 
quantum interference effects are investigated. For classical light, the photonic system realizes an 
optical analogue of coherent population trapping in the continuum encountered in atomic physics, 
where destructive interference between different light leakage channels leads to the appearance of 
a trapped state embedded in the continuum. For nonclassical light, two-photon interference effects 
are predicted, such as the tendency of photon pairs to bunch when decaying into the continuum. 
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I. INTRODUCTION 



It is well established that quantum mechanics of a non- 
relativistic particle and paraxial wave optics in dielectric 
media shear strong formal similarities (see, for instance, 
[ll-Q). Owing to such similarities, the temporal evolution 
of a quantum particle (e.g. an electron in an atom or in 
a crystal) can be mimicked by means of light propaga- 
tion in suitably-designed photonic structures. Quantum- 
optical analogies have seen in recent years a renewed 
and increasing interest, both theoretically and experi- 
mentally, mainly motivated by the possibility offered by 
optics to visualize at a macroscopic level certain coherent 
phenomena, originally proposed for quantum systems, 
which may be of difficult access or of controversial in- 
terpretation in the quantum context (see, e.g., [HHI and 
references therein). In particular, recent works theoreti- 
cally proposed and experimentally demonstrated the op- 
tical analogues of some important coherent phenomena 
encountered in atomic and molecular physics, such as 
coherent population transfer 0, 0| , electromagnetically- 
induced trans par ency [§| , and population trapping in the 
continuum [91-Hlj. The latter effect has been extensively 
investigated in atomic physics in connection, for instance, 
with the problem of autoionization of an atom by ultra- 
violet radiation or in laser-induced continuum structures 
(see [H], HH an d references therein) . In the simplest case, 
population trapping in the continuum describes decay 
suppression of two discrete states coupled to a common 
and unstructured continuum: as a single bound state 
coupled to the continuum decays in an irreversible way 
and population is transferred into the continuum, un- 
der certain conditions coupling of the continuum with 
the additional other bound state may partially or totally 
suppress the decay of states owing to a destructive inter- 
ference effect which was first studied by Fano for radiative 
transitions to autoionizing states in atoms [l4{ . Fano-like 
resonances in photonic systems have been considered in 
several recent works as well [15l 1 1 61 ] , with a main focus 
on the scattering (transmission) properties of the struc- 
tures. These previous studies in optics have mostly con- 



sidered propagation of classical light in coupled guiding 
structures, disregarding the quantum nature of light. For 
classical light, destructive interference among different 
tunneling paths is responsible for the existence of bound 
states in the continuum, similarly to Fano interference in 
atomic physics. However, when few photons or nonclas- 
sical beams are used to excite the photonic structures, 
light propagation may show interference effects of quan- 
tum nature. Since long time, coupled waveguides are 
known to behave similarly to beam splitters (see, for in- 
stance, [17]) and to show strictly quantum features when 
nonclassical light propagates through them, such as two- 
photon Hong-Ou-Mandel quantum interference originally 
demonstrated for beam splitters [18| . Like beam splitters 
[19j |. coupled waveguides offer the possibility to transfer 
nonclassical properties of light and to generate entan- 
gled states [171 [20] • With recent technological advances 
in the realization of high-quality low-loss integrated pho- 
tonic structures and nonclassical light sources, such pos- 
sibilities are nowadays realities. High-fidelity Hong-Ou- 
Mandel quantum interference and integrated optical re- 
alizations of other key quantum photonic circuits have 
been recently demonstrated in silicon-based waveguide 
couplers [2l| . In another recent experiment, Bromberg et 
al. |22| showed nontrivial photon correlations in coupled 
waveguide arrays and observed them using classical in- 
tensity correlation measurements. The possibility offered 
by experimentalists to test quantum aspects of light in 
integrated optic networks motivates to extend previous 
studies of quantum-optical analogies from the classical 
to the quantum level, highlighting interference effects of 
purely quantum nature. For instance, a recent theoretical 
study of optical Bloch oscillations and Zcncr tunneling in 
optical lattices [23[ showed that propagation of nonclas- 
sical light in the lattice may reveal certain particle-like 
aspects of light and quantum interference phenomena. 
In this work we investigate classical and quantum inter- 
ference effects of light in an optical analogue of coher- 
ent population trapping, based on two optical channel 
waveguides side coupled to a common slab waveguide re- 
cently proposed in Ref.[9j. For classical light (coherent 
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states), such system realizes in optics a simple analogue 
of destructive Fano interference with suppression of light 
leakage from the channel waveguides (analogous to two 
discrete states) into the slab waveguide (the analogue of 
the continuum) . When the description of light propaga- 
tion is made at the quantum level, photon pair excitation 
reveals strictly quantum features, such as the tendency 
of photons to bunch when decaying into the continuum. 



II. THE OPTICAL ANALOGUE OF 
POPULATION TRAPPING IN THE 
CONTINUUM: BASIC MODEL AND CLASSIC 
WAVE OPTICS DESCRIPTION 

A. The model 

Let us consider propagation of quasi-monochromatic 
and paraxial light wave packets at carrier frequency 
uj = 2ttc I A in a weakling guiding dielectric structure, 
with optical axis z and refractive index profile n(x,y), 
composed by two straight and parallel single-mode and 
equal channel waveguides Wi and W2, side-coupled to 
a slab waveguide S as shown in Fig. 1(a). In the parax- 
ial, weak guidance and quasi-monochromatic approxima- 
tions, the electric field can be written as E(x,y,z,t) oc 
ip{x, y, z,t) exp(j/3z — iut) + ip ji (x, y,z,t) exp(—i(3z + iujt), 
where the envelope tp varies slowly with respect to z and 
t over one wavelength A and one optical cycle 2tt/uj. Ne- 
glecting nonlinearities and group- velocity dispersion, the 
slow evolution of the electric field envelope ip(x,y,z,t) 
along the paraxial z direction is governed by the scalar 
wave equation 



i I ij) z + —ip t 



■■^n(lpxx+i>yy)+V(x,y)<lp, (1) 



where (3 — (u)/co)n s is the reference propagation 
constant, n s is the substrate refractive index, v g — 
(dft/duj)^ 1 is the group velocity of light, and V(x,y) — 
P[n s — n(x, y)\/n s . The electric field envelope ip has been 
normalized such that the cycle-averaged total energy of 
the electromagnetic field (see, e.g., [24j) is given by 



U 



dxdydzij^ 'ip — hw dxdydtij^ ip . (2) 



At the entrance plane z = 0, light is typically injected 
into either one, or in both, the channel waveguides Wi 
and W2, in the form of either monochromatic waves 
or wave packets with spatial profiles tailored to match 
their fundamental modes. Owing to evanescent coupling 
with the slab waveguide S, light leakage into the contin- 
uous set of modes of the slab is generally found, how- 
ever for certain geometric settings a trapping state may 
exist, corresponding to destructive interference between 
different tunneling paths into the continuum Q. At 
a classical level, light propagation in the waveguides is 



well captured in the framework of a coupled mode equa- 
tion approach, as previously shown in Rcf49]. After in- 
troduction of the spectral decomposition ip(x,y,z,t) = 
(27T)- 1 / 2 / dQct>(x,y,z,tt)exp(-int), from Eq.(l) it fol- 
lows that the spectral amplitude 4>(x, y, z, tt) satisfies the 
wave equation 

^* = -^ai&xx + <i>yy) + V(x,y)(j>- —(f). (3) 

Let us indicate by U\(p) and u 2 (p) the fundamental 
modes of waveguides Wi and W2 , and by u k (p) the con- 
tinuous set of modes of the slab waveguide S with the 
normalization conditions J dp\u±(p)\ 2 — J dp\u2{p)\ 2 — 
1, / dpuk{p)u* k ,{p) = 5{k — ft') and / dpui{p)u* 2 {p) ^ 
J dpui(p)u* k (p) ~ / dpu 2 (p)u* k (p) ~ 0, where p = (x,y). 
Let us search for a solution to Eq.(3) in the form 



<j>(p,z,Q) = [ct(z, n)tti(/?) + c 2 (z,Cl)u 2 (p)- 



(4) 



+ 



J dkc k (z,Q)u k (p) 



exp(i£lz/v g ). 



The evolution equations of modal amplitudes c\ , C2 and 
Cfc then read Q 



.dci 

' dz 



A/3 ci + J dkgi(k)c k 

dc f 
i— = A(3 c 2 + dkg 2 (k)c k 



(5) 
(6) 



= A/3(k)ck + 9*(k)ci + g 2 (k)c2 (7) 

where: gi(k) is the coupling amplitude between modes U\ 
and u k ; g 2 (k) is the coupling amplitude between modes 
u 2 and Uk] and A/?(fc), A/3o are the propagation constant 
shifts of modes u k and ui <2 , respectively, from the refer- 
ence value (3. In their present form, Eqs.(5-7) are analo- 
gous to the dynamical equations describing the quantum 
mechanical decay of two bound states Wi and W2 into 
a common continuum S [see Fig. 1(b)], provided that the 
temporal dynamics of the quantum mechanical problem 
is replaced by the paraxial propagation in space of light 
waves. The two bound states, with the same energy, are 
embedded in the continuum provided that A/3o falls in- 
side the continuous spectrum Af3(k), a condition which 
is satisfied whenever the refractive index change Arts in 
the slab waveguide is smaller than the index change An g 
in the two channel waveguides (see @ for more details). 
Note that, if waveguides Wi and W2 are symmetrically 
placed at opposite sides from the slab waveguide S as 
shown in Fig. 1(a), the coupling coefficients gi(k) and 
g 2 (k) are the same, i.e. g 2 (k) = gi(k) = g(k), and in this 
case a trapped state does exist, as discussed in Ref.[§]. 
Although trapped states in the continuum may also ex- 
ist when the waveguides Wi and W2 are horizontally 
displaced or placed at different distances from the slab 
waveguide Q , in this work we will consider the simplest 
symmetric case shown in Fig. 1(a). 
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B. Optical analogue of population trapping 

If one of the two channel waveguides, say e.g. W2, 
were removed, light initially injected into waveguide Wi 
would decay into the slab waveguide, a phenomenon fully 
analogous to the quantum mechanical decay of a bound 
state coupled to a continuum. In the markovian approxi- 
mation, valid for weak coupling and for a nearly unstruc- 
tured continuum, the decay is well described by an ex- 
ponential law. The presence of waveguide W2 generally 
modifies the decay behavior and, under certain condi- 
tions, the decay can be suppressed owing to a destructive 
Fano-like interference of different decay channels. Full or 
fractional suppression of the decay is related to the ap- 
pearance of a trapped (or dark) state in the continuum. 
To derive the decay laws of light waves in the channel 
waveguides Wi and W2 , we follow a standard procedure 
[HI, detailed for instance in Ref.@, and eliminate the 
amplitudes Ck of continuous modes by a formal integra- 
tion of Eq.(7) with the initial condition Cfe(0,fi) = 0. 
This yields a set of two coupled integro-differential equa- 
tions for the amplitudes c± and C2 of discrete modes. In 
the weak coupling limit and assuming a nearly unstruc- 
tured continuum, the markovian approximation can be 
made and the following reduced equations for amplitudes 
C\ and C2 are derived 

3c\ dc 2 

— = -~iA/3 c 1 -a(c 1 +c 2 ) , — = -iAf3 c 2 -<7(ci+c 2 ), 

(8) 

where 

a = L dr / dfc lf( fc )| 2ex P{-*[ A ^( fc ) - A A)]t} (9) 

is the decay rate of the single channel waveguide into the 
continuum. The solution to Eqs.(8) reads explicitly 

ci(*,n) = SnOOci(o,n) + Si 2 (*)c2(o,n) (10) 

c 2 (z,il) = S 21 (z) Cl (0,Q.) + S 22 (z)c 2 (0,O) (11) 

where 



Sn(z) = - cxp(-iAf3 z) [1 + cx P (-2ctz)] (12) 



■ exp(-iA(3 z) [1 - exp(-2erz)] (13) 



S 22 (z) = Sn(z) , S 2 i(z) = Si 2 (z). 



(14) 



Note that, in the quasi-monochromatic approximation 
assumed in this work, the matrix coefficients S nj i are in- 
dependent of frequency Q. To understand the appearance 



of the optical analogue of population trapping, let us con- 
sider the monochromatic case, with the only nonvanish- 
ing spectral component at ft = 0, and two different input 
excitations, corresponding the former to single waveguide 
excitation and the latter to simultaneous excitation of 
the two channel waveguides. In the former case, assum- 
ing for instance ci(0,Q) = \/2tt6(Q) and c 2 (0,Q) = 0, 
one obtains 



rl>(p,z) = S 11 {z)u 1 (p) + S 12 (z)u 2 {p) + S 13 (z)0(p,z) (15) 
where the last term on the right hand side of Eq.(15) 
accounts for the light field tunnelled into the slab waveg- 
uide and 9{p, z) ia a suitable superposition of continuous 
modes Uk(p) normalized such that J dp\9(p, z)\ 2 — 1. For 
power conservation, the relation |5n | 2 + IS12I 2 + |<Si3| 2 = 
1 then holds. Note that, after a propagation distance 
z a few times the decay length Id = (1/er), one has 
\Sn\ 2 = \S 12 \ 2 = 1/4 and |S 13 | 2 ~ 1/2, i.e. half 
of the injected light power has decayed into the slab, 
whereas the other half of light power is equally dis- 
tributed into the two channel waveguides. The fact that 
the decay is not complete (fractional decay) indicates 
that a bound state embedded in the continuum does ex- 
ist. When both waveguides Wi and W2 are initially ex- 
cited with coherent fields of equal amplitudes but op- 
posite sign, i.e. Ci(0,Q) = — c 2 (0,Cl) — V2nS(fl), one 
obtains ip(Pi z ) — l u i(p) ~ u 2{p)] exp(— iAfi^z), i.e. the 
decay into the slab waveguide is fully suppressed. This 
is due to a destructive interference effect between dif- 
ferent decay channels when C2 = —C\ [see Eqs.(7) and 
(8)] and to the existence of a trapped state embedded in 
the continuum. Numerical examples of fractional light 
decay for single waveguide excitation, and of full decay 
suppression for simultaneous waveguide excitation in the 
trapped state, as obtained by a direct numerical analy- 
sis of Eq.(l) in the monochromatic regime, are shown in 
Fig. 2. In the simulations, we assumed circular channel 
waveguides with a Gaussian index core profile of radius 
r c (at 1/e), and a step-index slab waveguide of thick- 
ness a. Equation (1) has been integrated by a standard 
split-step pseudospectral method with absorbing bound- 
ary conditions [9j. 

Generalization of light propagation in the non- 
monochromatic case simply follows from the superposi- 
tion principle. For instance, if waveguides Wi and W2 
are excited at the input plane by two pulses with en- 
velopes ri(t) and r 2 (t), from Eqs.(4), (10) and (11) it 
follows that the field envelope ip(p,z,t) can be cast in 
the form 
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ip(p,z,i) 



Sn(z)n [t- — ) + S 21 {z)r 2 it- — 



Sii{z)n ( t - — ) + S 22 (z)r 2 ( t 



u 2 (p) + i>s 



(16) 



where ips — ips{p,z,t) is the fractional part of the field 
tunnelled into the slab waveguide. In particular, let 
us consider two coherent pulses with the same envelope 
but phase reversed and delayed by an interval 6, i.e. 
r\(t) = r(t) and r 2 (t) = —r(t — S). In this case, in- 
terference effects, leading to light trapping in waveguides 
Wi and W2 and full suppression of leakage into the slab 
waveguide S (ips — 0), occurs for 5 = or, approxi- 
mately, for a delay 6 much smaller than the characteris- 
tic pulse duration r c . For S larger than r c , the two wave 
packets are not overlapped and behave as independent 
beams, leading to fractional decay. The latter scenario is 
also observed if the two pulses are temporally overlapped 
(5 = 0) but incoherent, i.e. their phase difference changes 
randomly in time. 



III. QUANTUM DESCRIPTION AND 
NONCLASSICAL EFFECTS 

A. Quantization procedure 

To describe propagation of nonclassical light in the 
coupled waveguide system, the classical paraxial wave 
field ip in Eq.(l) or, similarly, the classical c-numbers Ci, 
c 2 and Ck in Eqs.(5-7), have to be replaced by quantum- 
mechanical operators satisfying suitable commutation re- 
lations, and different quantization procedures may be 
adopted. A first approach, which is well suited when the 
classical problem is formulated in terms of coupled-mode 
equations (5-7), is the input-output operator formalism 
commonly used for linear quantum-optical networks, ei- 
ther in the Heisenberg or in the Schrddinger pictures (see, 
for instance, (T3, [Hi). A second approach, suited when 
the classical problem in Eq.(l) is formulated as a prop- 
agative (rather than as an initial-value) problem, is to 
adopt a quantization procedure for the classical field ip 
as an evolution in space (rather than in time). In the 
Schrodinger picture, this leads to an evolution in space 
of a many-photon probability amplitude. Such a phe- 
nomenological approach, which will be adopted in the 
following analysis, has received a growing use and appre- 
ciation in quantum theories of optical solitons [26j-|28| ; its 
consistency with standard canonical quantization proce- 
dure has been discussed in Rcf.^. Similar procedures 
have been also developed to study, in the Heisenberg pic- 
ture, paraxial propagation of nonclassical light and ap- 
plied to problems of quantum imaging (see, for instance, 
[29| ) . The quantization procedure consists in writing the 



classical paraxial wave equation (1) in Hamiltonian form 
assuming the paraxial spatial coordinate z as an indepen- 
dent variable [2t| . Introducing the new field II = ihip^ 
and the Hamiltonian H = J dxdydtH with density 

H = -^{U x ip x + Uyi/ty) - — Ityt - iV(x, y)Thl>, (17) 

it readily follows that the Hamilton equations ip z = 
(SH/5U), U z = -(SH/Stp) yield Eq.(l) and its com- 
plex conjugate, so that n is canonically conjugated 
to ip. Quantization is then accomplished by replac- 
ing the classical fields ip and n with the operators 
ip(x,y,t) and II = ihip^ (x,y,t) satisfying the commu- 
tation relations [ip(p, t), ip^(p', t')] = 5(p - p')8{t - t') 
and [i>(p,t),i>(p',t')} = [f(p,t),f(p',t')} = 0, where we 
have set p = (x, y) . By introducing the spectral decom- 
position ip(p,t) — (2tt)^ 1 / 2 J cttl0(p, ft) exp(— iClt), the 
second-quantized Hamiltonian operator reads 



H = h dpdVL 



2/3 1 



+ ■ 



'yVyJ 



V{ P ) 





Vn 



Note that in the spectral domain the following commuta- 
tion relations hold for the operators (p(p, fl) and <jy(p, £1) 



fa, n),tf( P ', n')J = 5( P - P ')S(n - «')> 
4>(p,n),4>(p',n')\ = \&(p,n),fr(j,tfj\ =o. (19) 



Note also that the field energy U [Eq.(2)] corresponds 
to the operator U — Hu> J dpdil^ (p,fl)4>(p,tt). In the 
Schrodinger picture, the quantum field is described by a 
vector state | Q(z)) which evolves according to 



d\0 

dz 



H\Q) 



(20) 



whereas the operators (p(p,Vl) do not evolve with z [301 ] . 
The state |Q) can be expanded in Fock space as \Q) = 
°ti|/ (<!) ^ z )); where the n-photon number state 
|/(™)) is defined by (see, for instance, [261 ] ) 



= / dqdtt 



,/ (n) (q.n,«) 



0t( Pl ,ni)....0t( ftl ,n n )|o> 

(21) 

and where we have set q = (p\,p2,...), fi = 
(Oi, •■•)■ The normalization conditions J2 n l a «| 2 = 1 
and / dqdfl\f^{q,n,z)\ 2 = 1 are also assumed, which 
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ensure that (Q\Q) = 1. Note that U\f {n) ) = nhiu\f^), 
i.e. the Fock state I/*-™-*) is obtained from the vacuum 
state |0) by creating n photons with space-frequency 
weighting function f^ n \ The evolution equation for 
the weighting function is obtained by substitut- 

ing Eqs.(18) and (21) into the Scrodinger equation (20) 
and using the commutation relations of field operators 
[Eqs.(19)]. One then obtains 

(22) 

Owing to the bosonic nature of photons, solely symmetric 
functions f( n > should be considered. 

B. Nonclassical effects with monochromatic beams 

Let us consider first the propagation of monochro- 
matic fields, so that in Eqs.(18), (19), (21) and (22) we 
may disregard integration over different spectral com- 
ponents f2 and use a single renormalized bosonic cre- 
ation operator 4>'(p) at frequency Q = satisfying the 
commutation relations [3l[ [$(p), ${p')\ — S(p — p') 
and [<j>{p), <t>(p')] — w{p)i<ft{p')] = , which replace 
Eq.(19). The simplest rt-photon number state, denoted 
by \g)n, is obtained by assuming in Eq.(22) /^(q, z) = 
g(pi, z)g(p2, z)...g(p n , z), where the function g{p,z) sat- 
isfies the classical wave optics equation 

igz = ~2p(9xx +9yy) + V(x, y)g (23) 

with the normalization J dxdy\g(x,y, z)\ 2 = 1. In this 
case one has 

\g) n = ^=(J dxdyg(x,y,z)^(x,y)\ |0). (24) 

Physically, this state describes the excitation of the op- 
tical system with an n-photon number state input beam 
with a spatial profile g(x,y,0) at the entrance plane 
z = 0. The classic wave optics description of light prop- 
agation, discussed in Sec. II, is attained by considering a 
superposition of photon number states \g n ) with a Pois- 
sonian distribution with c-number a, i.e. the coherent 
state |(?;a) co h defined by 

\g;a) coh = ± expH j 2/2K (/ dp 9 ( P ,z^(p)) n 10). 

n— 

(25) 

One can readily show that the coherent state \g; a) co h is 
an eigenstate of the field annihilation operator (j>(x, y) 
with eigenvalue ag(x,y,z), i.e. (j>(x, y)\g; a) cob = 
ag(x,y, z)\g;a)coh, where g{x,y,z) evolves according to 
the classical wave equation (23). Therefore the expecta- 
tion value coh (g; a\cf>(x, y)\g; a) coh = ag(x, y, z) yields the 
classical solution of the wave equation (23) for an input 



beam profile ag(x, y, 0). More generally, for a nonclassi- 
cal state \Q) obtained by an arbitrary superposition of 
photon number states \g) n with amplitudes a n , one can 
readily show that the expectation value of (jy(x, y)(j)(x, y) 
yields the classic wave optics intensity distribution, 
namely {Q\<ft{x,y)<j>(x,y)\Q) = (n)\g(x,y, z)\ 2 , where 
(n) = X)n n l a «| 2 i s the mean photon number of the in- 
put beam. The quantum aspects of nonclassical light for 
single beam excitation may be revealed when the statis- 
tics of photons trapped in waveguides Wi and W2 are 
considered. For example, let us consider excitation of 
waveguide Wi in its fundamental mode at the input plane 
z — 0, so that g(x,y, 0) = Ui(x,y), and let us compare 
the statistics of photons that remain in waveguide Wi 
when the input beam is a photon number state (nonclas- 
sical light) or a coherent state (classical light). According 
to the analysis of Sec. II. B, the wave amplitude g evolves 
according to [see Eq.(15)] 

g(p, z) = Sn(z)ui(p) + Suu 2 (p) + Si 3 {z)0(p, z), (26) 

where l^nl 2 + \S 12 \ 2 + \S 13 \ 2 = 1 and 9{p, z) defines the 
normalized spatial profile of the field tunnelled into the 
slab waveguide. If we introduce the operators 

at . (27) 

4 = / ' dxdyu 2 {x 1 y)4> i {x,y) (28) 

4 = J dxdy9(x,y,z)4 l ' : '(x,y) (29) 

the commutation relations [oi,ot] = <5i,fe and [tt|,ot] = 
[di, dk] = (i, k = 1, 2, 3) hold. Assuming that the vector 
state \Q) is given by a superposition of photon number 
states \g) n with amplitudes a n , one can write 

00 

1 2) = E T^t ( 5llflt i + Sl2a2 + Sl3a *) " |0) ' (30) 

The joint photon distribution P{n\; n^, z) to find ri\ 
photons in waveguide Wi, n 2 photons in waveguide W2 
and n 3 photons in the slab waveguide S is given by 

P(n 1 ,n 2 ,n 3 ;z) = \{n u n 2 ,n 3 \Q)\ 2 (31) 

where we have set 

\ ni ,n 2 ,n 3 ) = 1 a^a^a^lO). (32) 
V ' n\\n 2 \n 3 \ 

The explicit expression of P(ni,n 2l n 3 ; z) can be ob- 
tained by a double binomial expansion of the opera- 
tor (S n a\ + S 12 a\ + S 13 al) n entering in Eq.(30). The 
marginal photon distribution Pi(ni,z) to find n\ pho- 
tons in waveguide Wi is then obtained as P\{n\\z) = 
Sr 2 ,n 3 =o- P ( n i' n 2,"3;^), and reads explicitly (see, for 
instance, [l7| ) 

Pi(m;z)= f] K| 2 ^J^(l-«) n - m (33) 
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where n = |S , n(z)| 2 . Similar expressions are obtained for 
the marginal photon distributions P2(^2, z ) and ^3(^3, z) 
for waveguides W2 and S after setting k — \ Si 2 {z)\ 2 and 
K = \Si3(z)\ 2 , respectively. Let us now suppose that the 
input beam is a coherent state with mean photon number 
(ri) = no and Poissonian distribution 



n 1 ^ exp(-n ) 



(34) 



From Eq.(33) it follows that the marginal photon distri- 
bution P\(ri\\ z) remains Poissonian with mean photon 
number (rii) = no|<Sii(.2)| 2 that decays along the prop- 
agation distance z according to the classical decay law 
given by Eq.(12). Conversely, if the input beam is a pho- 



ton number state, i.e. a r , 



from Eq.(33) it follows 



that the marginal photon distribution Pi(ni,z) is given 
by the binomial distribution 

P 1 (n 1 ;z) =(11) " (ni ~ Uo) (35) 

[Pi(ni; z) = for n\ > uq] with photon mean 
(rii) = n |S'ii(z)| 2 and variance (An 2 ) = |S , n(z)| 2 [l - 



|S'ii(z)| 2 ]rio- The binomial photon distribution highlights 
the particle-like behavior of photons undergoing the de- 
cay from waveguide Wi and it is analogous to that cre- 
ated by a beam splitter when excited by a photon num- 
ber state in one port, and the vacuum state in the other 
one (see, for instance, [H|). The photons in the initially 
excited waveguide behave like independent classical par- 
ticles and tunnel into the other waveguides is ruled by a 
Bernoulli trial (a coin toss) with a cumulative tunneling 
probability given by 1 — |Sii(z)| 2 . 

When the two waveguides Wi and W2 are simultaneously 
excited by two nonclassical independent beams, quan- 
tum signatures of light propagation resulting from quan- 
tum interference can be detected by photon coincidence 
measurements. To describe propagation of independent 
beams, let us notice that for a given set of normalized 
and orthogonal solutions gi(p, z), g 2 (p, z), ... to the clas- 
sic wave equation (23), one can construct the n-photon 
number state 



0ii!n2!. 



dpgi{p,z)ft(p) ] x [ dpg 2 (p, z)ft (p)) x ...|0) 



(36) 



which describes excitation of the optical system with a 
set of independent beams with spatial profiles 171, g 2 
carrying n\, U2,.-. photons (n — ni + n 2 + ■•■)• In par- 
ticular, let us assume that at the input plane the waveg- 
uides Wi and W2 are excited in their fundamental modes 
by single photon number states, i.e. gi(p, 0) = u\{p), 
52(10,0) = u 2 (p) and 



\Q(z = o)) = |ui,u2>i,i = aJa£|o), 



(37) 



where the operators a\ and a\ are defined by Eqs.(27) 
and (28). According to the classical analysis of Sec. II. B, 
the waves g\ and g 2 evolve according to 

gi(p,z) = Sn(z)ui(p) + Si 2 (z)u 2 (p) + S 13 (z)9(p, z) 

(38) 

g 2 (p,z) = S 2 i (z)ui (p) + S 22 (z)u 2 (p) + S 2 3(z)9(p,z) 

(39) 

with 523 = for symmetry reasons (32|, Si 2 = S21, 
Su = S 22 [seeEqs.(12-14)], and |5 n | 2 + 15 12 | 2 + 15 13 | 2 = 
1 for power conservation. From Eqs.(36), (38) and (39) it 
follows that the state vector of the system at the generic 



propagation distance z is given by 

\Q) = (Sua\ + S 12 al + Si 3 a3jJ (^S 2 ia\ + S 22 al + 5 23 a 3 

(40) 

where a 3 is defined in Eq.(29). It can be readily shown 
that the expectation value of the field intensity for the 
two-photon state (40) is given by the incoherent super- 
position 

(Q\$Hp)fcp)\Q) = \ 9 i(p,z)f + \g 2 (p,z)\ 2 , (41) 

the absence of interference being due to the lack of a defi- 
nite phase relationship between the two photons (see, for 
instance, (33j). The two independent beams which excite 
the two waveguides Wi and W2 behave, therefore, as two 
incoherent classical fields, and therefore according to the 
analysis of Sec. II. B half of the light power tunnels, on 
average, into the slab waveguide S for propagation dis- 
tances z 3> Id ■ As the operators a\, a\ and S3 create 
photons in waveguides Wi, W2 and W 3 , respectively [see 
Eqs. (27-29)], from Eq.(40) the joint photon distribution 
P(n\, n 2 ,ri3; z) can be readily calculated, and its nonva- 
nishing terms are explicitly given by 
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P(2,0,0;z) 
P(l,l,0;z) 
P(l,0,l;z) 
P(0,0,2;z) 



P(0,2,0;z) = 2\S n \ 2 \S 21 \ 2 = 1 [1 - exp(-4 ( rz)] 2 
\SnS 22 + S 12 S 21 \ 2 = 1 [1 + exp(-4 ( rz)] 2 

P(0, 1, 1; z) = |5u523 + 5 13 5 21 | 2 = i exp(-4 ( rz) [1 - exp(-4<7z)] 



2|5i 



31 I ^23 1 



X - [1 - exp(-4 ( 7z)] 2 



(42) 
(43) 
(44) 
(45) 



An important result is that, for propagation distances z 
larger than the characteristic decay length l d = (1/tr), 
the joint probability P(1,0, 1) to find one photon in 
waveguide Wx and one in the slab waveguide S vanishes, 
and similarly for P(0, 1,1) [see Eq.(44)]. Such a result 
indicates that, if one of the two initially injected photons 
tunnels into the slab waveguide, the other photon does 
the same, i.e. photon pairs bunch when decaying into 
the continuum. This is obviously a nonclassical effect 
similar to the two-photon Hong-Ou-Mandel quantum 
interference in a 50% beam splitter [3. In our case, 
the vanishing of P(1,0,1) [and similarly of P(0,1,1)] 
is related to a destructive interference between the 
probability amplitudes SuS 2 3 and S13S21 entering in 
Eq. (44) which describe two possible paths for the photon 
pair. 



C. Nonclassical effects with wave packets 

In the previous subsection, we have considered propa- 
gation of nonclassical light in the monochromatic limit, 
however in practice nonclassical light such as that gen- 
erated by spontaneous parametric down conversion in a 
nonlinear crystal is always polychromatic to some extent. 
Moreover, in experimental settings quantum interference 



such as the two-photon Hong-Ou-Mandel interference is 
typically revealed as a dip in intensity correlation (photon 
coincidence) measurements when the time delay between 
the two incoming wave packets is varied [l8|, [2l| . It is 
therefore useful to extend the analysis of Sec. III. B to the 
case of polychromatic wave packets. 
For the sake of definiteness, let us assume that waveg- 
uides Wi and W2 are excited in their fundamental spa- 
tial modes by two spectrally-narrow wave packets gen- 
erated by spontaneous parametric down conversion in 
type-I nonlinear crystal. In this case, the vector state 
\Q) at the input plane of the waveguiding system is a 
generalization of Eq.(37) and given by a superposition 
of two-photon states with a spectral function C(fii,f2 2 ) 
(see, for instance, 0, H3[) 

\Q(z = 0)) = J rffiidO 2 C(n 1 ,fi 2 )a^Oi)a 2 (fi 2 )|0) 

-t -t (46) 

where the operators d|(f2i) and d 2 (f2i) are defined by 

a{ 2 (n) = J dpu h2 (p)ft(p,n). (47) 

The evolution of vector state along the propagation dis- 
tance z can be obtained by a straightforward extension 
of the analysis of Sec. III. B and reads explicitly 



dSlidShCinu Sh) exp [i(fii + n 2 )z/v g ] S n (z)al(Qx) + S'i 2 (z)4(^i) + Si3(z)a 3 ( i) 
S 2 i(z)a\(n 2 ) + S 22 (z)a\(Q 2 ) + S 23 (z)al(n 2 )} |0) 



(48) 



where 

al(o) = J dpd(p,z)4>H P ,n). (49) 

To study quantum interference effects with wave pack- 
ets, and in particular the vanishing of the joint photon 



probability P(l,0,l;z) = P(0, 1,1; z) for z > l d found 
for monochromatic beams [Eq.(44)], let us introduce the 
integrated intensity correlation function at times t\ and 
t 2 in waveguides Wi and S defined as 



8 



P(h,t 2 ;z) = I d Pl dp 2 (Q\^{p 2 ,t 2 )^(p 1 ,t 1 )iP(p 1 ,t 1 )^( P2 ,t 2 )\Q), (50) 

J Ai J Ai 



where A\ and A 2 are two areas in the transverse (x, y) 
plane surrounding waveguides Wi and S, respectively. 
The correlation function V(ti,t 2 ) is proportional to the 
joint probability of detecting one photon in waveguide 
Wi at time t\, and one photon in waveguide S at time t 2 , 
after a propagation distance z from the input plane. In 
fact, (Q\$ (p 2 ,t 2 )fo (pi,t 1 )i)(p ll t 1 )iJ;(p 2 ,t 2 )\Q) is pro- 
portional to the joint probability of detecting one photon 
at point pi and time t\, and one photon at point p 2 and 
time t 2 at the same propagation distance z. The integral 
over the areas A\ and A 2 thus gives the joint probabil- 
ity of detecting one photon in waveguide Wi at time ti, 
and one photon in waveguide S at time t 2 . For photons 
generated by spontaneous parametric down conversion 
in type-I nonlinear crystal using a monochromatic pump 
beam at frequency uj p = 2u, the two-photon spectrum 
C(fli,Q 2 ) of a pair of signal and idler photons can be 
expressed as [3 H3] 

C(fi 1) 2 ) = 5(Oi + n 2 )C?(Oi ) 2 )exp 

(51) 



where r = t 2 —t\ and where we introduced the real-valued 
correlation function r(r) defined by 

r(r) = — / dQG(Q, -tt) exp(-iQr). (53) 
2n J 

In practice, coincidence measurements correspond to an 
integration of V(t±, t 2 ; z) with respect to the time dif- 



where a — drr 2 (r) and the expression of the co- 
efficients S\i{z) and Si 2 (z) are given by Eqs.(12) and 
(13). The behavior of V(S;z) versus S shows a charac- 



I 

where G(p,i,Cl 2 ) is the phase matching function. The 
last exponential term on the right hand side in Eq.(51) 
has been introduced to account for a possible time de- 
lay 6 between signal and idler wave packets introduced 
by different optical paths from the crystal to waveguides 
Wi and W 2 . The phase matching function G(£li,fl 2 ) 
is assumed to be a real- valued and symmetric function 
[i.e. G(rti,fl 2 ) = G(fi 2 ,Qi)], peaked at around Oi = 
f?2 = 0, with e.g. a Gaussian profile [l~8j]. Taking into 
account that i>(p,t) — (27r) -1 / 2 J dVt(f)(p, VL) exp(— iflt), 
substitution of Eqs.(48) and (51) into Eq.(50), using the 
commutaion relations [4>(p, CI), a\ 2 (S^')] = ui^ 2 (p)8(Cl — 

Of), [$(p,Q),al(ft')] = d(p,z)S{n - ft') and the re- 
lations J Al dp\ Ul (p)\ 2 ~ 1, J A2 dp\9(p,z)\ 2 ~ 1, 
J Ai2 dp\u 2 Ap)\ 2 - f Ai dp\6(p,z)\ 2 = J A2 dp\u 2 (p)\ 2 ~ 
0, after some lengthy but straightforward calculations 
one obtains 



(52) 



ference r = t 2 — t\ over the resolving coincidence time, 
which is typically much longer than the correlation time 
t c of <?(t). Integrating Eq.(52) with respect to r from 
—oo to oo and taking into account that S 2 \ = S\ 2 and 
|S 23 | 2 = |5i 3 | 2 = l-|5n| 2 -|5 12 | 2 , the following expres- 
sion for the correlation function V versus time delay 5 is 
finally obtained 



(54) 

I 

tcristic dip at S — of width ~ r c , which is analogous to 
the Hong-Ou-Mandcl dip observed in two-photon inter- 
ference from a beam splitter [1S|] . Far from the dip, 'P(S) 



J 



V(h,t 2 ;z) = \r(r + 5)S n {z)S 23 (z) + r(r - S)S 13 (z)S 21 (z)\ 2 

I 



I 

V(S; z) = a [1 - |S u (z)| 2 - \S 12 (z)\ 2 ] [|5 n (z)| 2 + \S 21 (z)\ 2 + 

f drrir — S)r(r + S) 
+ 2Re[S 11 (z)S* 21 (z)} J -°° ' - ' \ 

Loo dTr ( T ) 



9 



reaches a constant value a(l — \Sn\ 2 — |Si2| 2 )(|Sii | 2 + 
|5 2 i| 2 ). For z larger than the characteristic decay length 
Id = l/u, the correlation V(S;z) vanishes at the dip 
(5 = 0, i.e. when signal and idler wave packets are tempo- 
rally overlapped, according to the analysis of Sec. III. B. 

IV. CONCLUSIONS 

In this work we have theoretically investigated propa- 
gation of classical and nonclassical light in a waveguide- 
based photonic structure that provides an optical ana- 
logue of population trapping in the continuum encoun- 
tered in atomic physics. For classical light waves, 
coupled-mode equation analysis, previously studied in 
Ref.[9(, shows that Fano interference between different 
light leakage channels is responsible for the appearance 
of a trapped state embedded in the continuum. To study 
propagation of nonclassical light, a second quantization 



model for the scalar wave equation, in the paraxial and 
quasi-monochromatic approximations, has been adopted. 
As for input beam excitation in a coherent state the clas- 
sical picture of light propagation is retrieved, quantum 
interference effects with no classical counterpart have 
been highlighted for photon number state excitation of 
the waveguide structure. In particular, the tendency 
of photon pairs to bunch when decaying into the con- 
tinuum has been predicted. Such an effect, which is 
similar to the two-photon Hong-Ou-Mandel interference 
in a beam splitter il.8|, may be observed as a dip in 
photon coincidence measurements. Our results indicate 
that photonic structures originally designed to mimic 
with optical waves the classical analogues of quantum- 
mechanical phenomena encountered in atomic, molecular 
or condensed-matter physics |5| , may exhibit themselves 
a strictly quantum behavior when single photon level is 
reached. 
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FIG. 1: (a) Schematic of the coupled waveguide system com- 
posed by two equal single-mode channel waveguides Wi and 
W2 side-coupled to a slab waveguide S. (b) The quantum 
mechanical analogue, showing the decay of two bound states 
coupled to a common continuum. 
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FIG. 2: (color online) Behavior of fractional light power 
trapped in channel waveguides Wi (blue curve) and W2 
(green curve) versus propagation distance in a L — 6-cm-long 
structure and corresponding transverse light intensity distri- 
butions at the input (z — 0) and output [z = 6 cm) planes, 
as obtained by numerical simulations of the wave equation 
(1) for two different excitation conditions. In (a) waveguide 
Wi is excited in its fundamental mode, leading to fractional 
decay. In (b) the trapped state is excited, corresponding to 
full suppression of light leakage into the slab S. Parameter 
values used in the simulations are: A = 980 nm, n s = 1.52, 



An s = 0.0049, 
a = 4 11m.. 
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